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Abstract 

We study a class of evolutionary partial differential systems with two components related to second order (in 
time) non-evolutionary equations of odd order in spatial variable. We develop the formal diagonalisation method 
in symbolic representation, which enables us to derive an explicit set of necessary conditions of existence of higher 
symmetries. Using these conditions we globally classify all such homogeneous integrable systems, i.e. systems which 
possess a hierarchy of infinitely many higher symmetries. 
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1 Introduction. 

Classification of integrable equations of a given family of nonlinear equations is an important topic in the field of 
soliton theory. There are many approaches to this problem, among which the symmetry approach has proved to 
be very efficient and powerful method. The symmetry approach has been used to classify large classes of integrable 
nonlinear partial differential equations and difference-differential equations such as scalar evolutionary equations, 
Volterra and Toda type equations, hyperbolic equations etc. We refer the reader to the recent review paper [T] and 
the references therein for details. 

This paper is devoted to a problem of classification of systems of integrable evolutionary equations with two compo- 
nents. The first work in this direction had been carried out by Mikhailov, Shabat and Yamilov They classified 
integrable second order systems of the form 

ut = A{u)u^^ + F{u,u^), det(yl(u)) ^ 0, u ^ (u, v f , 

where the 2x2 matrix A{u.) can be reduced to 

^(u) — diag(a, b), a = —1, 6=1, 

if the system possesses higher order conservation laws. Later, Svinolupov [5] studied the Burgers type equations, i.e. 
the case a = b = 1 (see also the review paper [B]). After these two special cases, attention was turned to the more 
general case: a ^ ±b and a6 ^ [7] and to the systems of higher order. In principle, the symmetry approach can be 
applied to classification of higher order systems. Although the method is algorithmic, the computation involved grows 
very rapidly with the order of the system, so further achievement in this direction had to wait for the development of 
the appropriate computer algebra. 

Recently, several research groups have developed the computer programmes to carry out the task of classification of 
polynomial integrable evolutionary systems with two components. Foursov and Giver considered the symmetrically- 
coupled evolutionary equations of lower order [5] . After linear transformation, the leading linear part of the systems 
can always be diagonalised. Tsuchida and Wolf studied polynomial homogeneous integrable evolutionary systems of 
mixed scalar and vector dependent variables [5] of order 2 and 3, where the authors compared their comprehensive 
results with the existing ones. The linear parts of systems they considered are diagonal including zero eigenvalues. 
Their classification methods have led to interesting new integrable equations, but they could not claim that their lists 
are complete in the sense that there are no other integrable systems of the fixed order in certain classes since the 
methods they used were based on searching symmetries of specific orders and forms. The symbolic representation 
combined with applications of number theory [TU] is the solution to this problem at least for subclass of homogeneous 
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polynomial evolutionary systems. It enables us to translate the question of integrability into the problem of divisibility 
of certain symmetric polynomials, which can be solved by using number theoretical methods. A series papers have 
been published on it including the complete classification of second order two components evolutionary systems 
and Bakirov-type systems [121 [13] . 

The ultimate goal of classification of integrable equations is to obtain the classification result for integrable equations 
of any order, i.e. the global classification result. So far the only global result has been obtained for scalar homogeneous 
evolutionary equations [141 115j . In this paper, based on the symbolic representation in |16| . we develop the formal 
diagonalisation method so that we can globally classify homogeneous two-component evolutionary systems whose 
leading linear terms are not explicitly diagonal. We demonstrate our method classifying integrable systems of the 
form 

\ = a^"+^"''u + F(u,Ma;,Ma:x,---,9^"-'~u,'(;,'ya;,'ya:x,---,9^""''w), n=l,2,3,..., T £ {0, 1 , . . . , n} , 

of arbitrary order (i.e. for any integer n > 0). Here is a homogeneous differential polynomial in u,v and their 
x-derivatives. If function F does not depend on v, Vx, ■ • ■ , d^T^v then we can exclude v from system ([ij and obtain 
non-evolutionary equations of odd order 

u« = 52"+iu + if(ii,w„u,,,...,af wt), K^DliF), 1,2,3,... (2) 

The even-order non-evolutionary equations of form ^ have been studied using the perturbative symmetry approach 
in symbolic representation in [16] . The famous Boussinesq equation [Til 

'^tt — ^xxxx ('^ ^xx 

belongs to this class. As a result, integrable equations of orders 4 and 6 were classified as well as three new integrable 
equations of order 10 were found. The approach developed in [16] is also suitable for odd order equations. However, 
there is a difficulty to formulate the explicit necessary conditions for this case. Besides, it can not be used for global 
classification. The special subclass of non-evolutionary equations was also considered in '20] , where a family of partially 
integrable equations with only one higher symmetry was constructed. 

In this paper in the framework of the perturbative symmetry approach in the symbolic representation we derive an 
explicit set of necessary conditions of existence of a hierarchy of higher symmetries for systems ([TJ. These conditions 
are obtained after diagonalisation of the linear part of system ([TJ via a non-local linear transformation of dependent 
variables u,w. We prove the ultimate classification result: If a homogeneous system ([iJ with positive weight of u 
possesses a hierarchy of infinitely many higher symmetries, then up to re-scaling it one of the following two equations: 

1 U4 = Uxx + iuVx -t- VUx — Su'^Ux, 

\ vt^ [dx + -v^, n = 1,2,3,.... 

These systems can be rewritten in the form of non-evolutionary equations if we introduce a new variable u = Wx- 

Wtt = Wxxx + iWxWt^x + WxxWt - iwlwxx- 

wtt = {dx + - wl 

The first equation is known to be integrable [TS]. The second equation can be brought into linear fu = 9^"^^/ by the 
Colc-Hopf transformation w — log(/). 

The paper is organized as follows: in section[2]we remind basic definitions and notations of the symbolic representation 
[16j : in section [3] we develop the formal diagonalisation formalism for systems ([1]); using the formal diagonalisation 
method we obtain explicit formulae for the (approximate) symmetries of system ([T]) and derive necessary conditions for 
the existence of an infinite hierarchy of higher symmetries; in section[2]we apply these conditions to classify integrable 
homogeneous systems ([T]) for all orders n = 1,2,3,... and for all possible representations r = 0, 1, . . . , n. 

2 Symbolic representation over the ring of differential polynomials. 

In this section we remind basic definitions and notations of the ring of differential polynomials and symbolic represen- 
tation (for more details see e.g. [HI [TBI [19] ) . 
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Throughout the whole paper we assume that all functions, such as F and K in equation ([T]) and are dif- 
ferential polynomials in variables u,v,Ux,Vx,Uxxt ■ ■ ■ ■ We introduce an infinite sequence of dynamical variables 
{uq, vo,ui,vi,U2,V2, ■ ■ ■ } by the identification 

Uq = U, Vq = V, Un = Vn = O^V. (3) 

We often omit the zero index and write u and v instead of uq and vq. 

We denote TZ the ring of polynomials over C of infinite number of dynamical variables. We also assume that 1 ^ TZ. 
Elements of the ring TZ are finite sums of monomials with complex coefficients and therefore each element depends on 
a finite number of the dynamical variables. The degree of a monomial is defined as a total power, i.e. the sum of all 
powers of dynamical variables that contribute to the monomial. Let TZ" denote the set of polynomials in TZ of degree 
n. The ring TZ has a natural gradation 

= ^" ' ^" • ^" C 7^"+'" . 

rieZ+ 

Elements of TZ^ are linear functions of the dynamical variables, TZ^ quadratic, etc. We suppose that 1^7?.. It is 
convenient to define a "little-oh" order symbol o{TZ"). We say that / = o(7^") if / £ ®fe>„''^'^j i-G- the degree of 
every monomial of / is bigger than n. 

Let /i, ly be two positive rational numbers0, which we call the weights of u and v respectively and denote W{u) = 
IJ,, W{v) = J/. We define the weights of dynamical variables ^ as W{ui) = n + i and W{vj) = v + j. The weight of a 
monomial is the sum of the weights of dynamical variables that contribute to the monomial. We say that a polynomial 
/ G 7^ is a homogeneous polynomial of weight A (and write W{f) = A) if every its monomial is of weight A. The 
weighted gradation of TZ is defined as 



TZ 7^m/^+ni/+s 1 TZp • TZq CZ TZp 



where TZp is the set of polynomials in TZ of weight p. In fact, we can define a degree- weighted gradation of the ring TZ 

n,m, sfEZ>o, rf^Z^- 

where TZp is the set of polynomials in TZ of degree r and weight p. Note that the condition > makes each such 
subspace to be finitely generated. For example, ii p — 1, v = 2 then 

TZl = span{uM2, uf, uiv, uvi, v^} 
The ring 7?, is a differential ring with a derivation defined as 

^^-Ef-^+ia^+-^-+i^)- (4) 



k>0 



Since 1^7^, the kernel of the linear map ■ TZ Im C 7?. is empty and therefore D,^, ^ is defined uniquely on 
Im Dx ■ It is easy to verify that 

To any element g ^ TZwe associate differential operators and g^^^ called Frechet derivatives with respect to u and 
V and defined as 

k>0 k>0 

Now we define the symbolic representation TZ of the ring TZ. A symbolic representation of a monomial 
Uq°u'1^ ■ ■ ■ UpP-u™"?;™! • • • w™' , no + ni H h rip = n, ttiq + nii -\ h = m 

is defined as: 

nn ni n„ mn mi rrin 

— > U W (Cl?2 ' ' ' + 1 ' ' ' Cio+rii ' ' ' ^n)? (Cl C2 ' ' ' CmoCjio + l ' ' ' Cmo+mi ' ' ' Cm)c ' (^) 



^In principle, the weights could be any rational numbers including zero, but in this paper we consider only positive weights. 
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where triangular brackets ()j and O^; denote the averaging over group I]„ of permutations of n elements ^i, . . . ,^„, , 
i.e. ^ 

(c(Cl, . . . ,Cn,Cl, • ■ ■,Cm))i = y2 C(cr(^i), . . . ,f7(C„),Cl, ■ ■ • , Cm) 

and group S„j of m elements Ci, ■ • ■ , Cm respectively. Later we refer to this as symmetrisation operation. For example, 
linear monomials u„ , are represented by 

Un ^ wCr, "m ^ ^^Cr (6) 

and quadratic monomials UnUm, UnVm^ VriVm have the fohowing symbols 

?/ 7; 
UnU^ -> y (CrCr + CrC2 ), «n^^m - "^^ICKD, «n^^m - y(Ci"C2"' + Q) ■ (7) 

To the sum of two elements of the ring corresponds the sum of their symbols. To the product of two elements f,g(zTZ 
with symbols / u"v"^a{^i, . . . , Ci, • • ■ , Cm) and g uPv'^b{^i, . . . ,^p,(i, ■ ■ ■ , Cg) corresponds: 

fg ^ «"+Pz;™+«((a(a, . . . , Cl, . . . , Cm)&(Cn+l, ■ • • , Cn+p, Cm+1, ■ ■ • , Cm+«))«)c, (8) 

where the symmetrisation operation is taken with respect to permutations of all arguments ^ and C- It is easy to see 
that the symbolic representations of quadratic ([7]) and general (O monomials immediately follows from ^ and 

li f G TZ has a symbol / u"u"'a(^i, . . . ,CniCii ■ • ■ 7 Cm), then the symbolic representation for its A^-th derivative 
D'^{f) is: 

D^{f) ^ + 6 + • • • + Cn + Cl + C2 + • ■ • Cm)'^a(Cl, . . . , Cn, Cl, • • ■ , Cm). 

We will assign a symbol 77 to the operator in the symbolic representation with obvious action rule 

77^(u"«™a(6, . . . , C„, Cl, • • • , Cm)) = ^/"z;"(Cl + 6 + • ■ • + Cn + Cl + C2 + • • ■ Cm)'^a(Cl, • • ■ , Cn , Cl , • • ■ , Cm) 

Note that N in the above formula can be not only positive integer but in principal any rational number (or even a 
complex number) as long as such calculation has any sense apart from formal . In the next section we will need N to 
be half-integer numbers which correspond to fractional derivation in a;-space. 

If g e 7^ and g —>■ M"u'"a„.„i(Ci, • . • , Cm Cij • ■ • 7 Cm) then for the symbol of its Frechet derivatives g^^u and 17,.^ we have 
nM""^i;'"a„,„(Ci, . . .,£,n-i,V,Ci, ■ ■ ■ ,Cm), g*,v ^ m-«"'(;™~^a„,m(Ci, • . • ,Cn,Ci, ■ • ■X^n-i,^)- 

Thus we obtained the symbolic representation TZ of the differential ring TZ. 



3 Structure of symmetries and approximate symmetries. 

Consider a family of homogeneous polynomial systems ([T]), i.e. 

Vt^dl"+'^^''u + F{u,Ui,...,U2n~r,V,Vi,...,V2„-r), n=l,2,3,..., r G {0, 1 , . . . , 7l} . 

Definition 1. A pair of differential polynomials G and Mis called a symmetry of an evolution system Ut = /, f t = g, 
where f, g, G and M are all in ring TZ, if system Ur — G, Vr — M is compatible with the given system. 

This definition is equivalent to the Lie bracket between a — (/, g)^"^ and b = (G, Af)^'' vanishing, where the bracket is 
defined as 

r , . UAG) + f*AM) ~ G.M) - G.Ag) \ 

^ ' ^ V 9*.u{G) + g*AM) - M^Af) - M.Ag) J ■ ^ ' 

The compatibility condition for system fO]) can be written as 

Dt{G) ^ D'^Vr ^ DIM, Dt(M)^Dr{u2n+l-r + F), (11) 

where evolutionary derivations Dt and Dj- are 

j>0 \ J J / j>0 \ J J 

Eliminating M from these equations, we obtain 

D-''Dj{G)=Dr{u2n+l-r+F), 

Thus the symmetry of system ^ is completely determined by its first component G G TZ. 
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Definition 2. We say that a differential polynomial G G 7^ generates an approximate symmetry of degree p of system 
(0j if G satisfies the equation: 

D-^DUG) - Dr{u2n+l-r + F) = o{nP)- 

Notice that any system ^ has an infinite hierarchy of approximate symmetries of degree 1. These are simply the 
symmetries of its hnear part Ut = Vr, Vt = U2n+i-r- The requirement of the existence of approximate symmetries of 
degree 2 is very restrictive. Integrable equations by definition possess infinite hierarchies of approximate symmetries 
of any degree. We wiU prove later in this paper that the existence of infinite hierarchy of approximate symmetries of 
degree 2 and the existence of at least one exact symmetry for systems ([SJ implies integrability just as in the case of 
scalar evolutionary equations [llj . 

We now derive the necessary and sufficient conditions of the existence of approximate symmetries of degree 2. It is 
convenient to do this in the symbolic representation. 

System ^ in the symbolic representation takes the form 

{ Vt^ <2„+l-. ^ j^^^^ j^k^^ u^v'^-^a,,k-,iCl, .■■,^^Xl,■■■Xk-^). ^^^^ 

Let us assume that the weight of variable u, denoted by w, is positive, i.e. W{u) = w > 0. Since the system is 
homogeneous, we have 

f w + W{Dt) = W{v) + r 

\ W{v) + W{Dt) = 2n + l- r + w, 

which leads to W{Dt) = n + | and W{v) = w + n + ^ ~ r. 

Homogeneous polynomials a^. j (^i , . . . , , Ci , . . . , ) are symmetric with respect to variables , . . . , and Ci , . . . , . 
The degrees of a^.j are given by 

2—7 

deg(aij) = 2n + l- r + u; - iW{u) - jW{v) = (2 - j)n ~ {i + j ~ l)w + {j - l)r H (13) 

The degree of aij determined above must be a non-negative integer, otherwise, we put aij = 0. The sum in (fT2|) 
terminates due to the assumption W{u) = w > 0. 

The symmetry of system ^ always starts with linear terms. A homogeneous symmetry of © in the symbolic 
representation starts either with m^™ or with uC™^''. Without loss of generality we have 

s 

s>2 j=0 

or 

s 

Ur^G^ + E E ^'v^'^A,,,.,{^,,. . . Ci, . . . , Cs-j), m > (15) 

s>2 j"=0 

We call the symmetries of the form (ITil) as type I symmetries and those of the form (fTS|) - type II symmetries. 
Without causing confusion, we call integer m in (|14p or psp the order of the corresponding symmetry. Functions 
Ai,j (^1 , . . . , ^i, Ci , . . . , Cj) in (fT4| and (fT5|) are homogeneous polynomials in their variables, symmetric with respect to 
arguments , . . . , and (^i, . . . ,Q- These functions can be explicitly determined in the terms of system ()12p from the 
compatibility conditions. 

3.1 Formal diagonalisation. 

Let us first concentrate on how to compute the Lie bracket defined as (fT(I| between the linear part of system (fT^ 
denoted by K^, i.e. 

k'.(^z^l:)^l(i). L.(^.z-r"i) (16) 

and any pair of differential polynomials. We know its symbolic representation takes simple and elegant form if matrix 
L is diagonal [TT]. Inspired by this, we shall formally diagonalise matrix L, produce the required formula in new 
variables and then transform back to the original variables. 
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Notice that matrix L has two eigenvalues ±77"+ 2 . Therefore, there exists a hnear transformation 



T 



such that 

Let us introduce new variables u and ij 



1 1 



T-^LT = diag(?7"+^ -77"+^) 



V \ V 



u + i) 



Equally, we have 



w y 2 V u-r;-"-2+'-(t;) 

The new variables w and v have the same weights, i.e. W{u) = = W{u). Without causing confusion we assign 

the same symbols ^ and C for the symbolic representation of the ring generated by u, v and their derivatives. The 
exponents of symbols can be half-integer, which corresponds to half-differentiation in x-space. 

Before we work out the exact form of system ()12p in variables ii and v, we prove two useful propositions. 

Proposition 1. Under the transformation T , an arbitrary polynomial . . . , ^i, C17 • ■ • 7 Cj) takes the form 

uV6,,,(a, . . . ,C„Ci, . . . ,0) = ^^K^+^f''+^-P-^CfCj(-l)^-« (17) 

(Cl I • • ■ I ?P7 Cl I ■ ■ • J Ci-pi ?p+l I • • ■ I ?p+gj Ci-p+1 1 ■ ■ • J Ci+J— p— g) (Cp+1 ' ' ■ Cp+ijCi— p+1 ' ■ ' Ci+j-p-<j) 

where Cl are binomial coefficients defined by Cf — jrp^jjT- 

Proof. We prove the statement by induction on both i and j. It is easy to see that linear terms u^l, v(f transform as 

The procedure is: we first substitute ^1 and Ci by 77, the symbol of the total x-derivative D^, u and v hy u + i) and 

u^^ ^ — v(^i ^ ; and then we compute the action of 77 on u -I- -0 and u£,-^ ^ — v(^-^ ^ . For the higher degree 
terms, we apply this procedure for every element in the arguments. Assume formula (|17p without symmetrisation is 
true for arbitrary i and j . We have 

p=0 q=0 

bi+lj{(,l, ■ ■ ■ ,S,p,ClT ■ ■ J Ci-pj ^p+1: ■ • ■ I Cp+gi Ci-p+li ■ • ■ 7 Ci+j~p-q){'^ + v) 

bi+lj {£,1 J ■ ■ ■ I Cpi Cl I ■ ■ • I Ci— p J ^p+1 J • ■ ■ I ^p+g+1 J Ci— p+1 : ■ ■ • I Ci+j—p — q) 
p=0 9=0 

bi+lj {£,1 J ■ ■ ■ I Cpi Cl 5 ■ ■ • 5 Ci— p J Ci+1— pi Cp+1 1 • ■ ■ J Cp+<3i Ci— p+2 j • ■ ■ j Ci+J— p— 9+1 ) 

= E E ■ ■ • ep+90-p+2 • • • C.+ l+,-p-9)"+^"' 

p=0 9=0 

bi+l.j (Cl 5 ■ ■ ■ I Cpi Cl : ■ ■ • I Ci+1— pi Cp+1 5 ■ ■ ■ I Cp+9i Ci— p+2 : ■ ■ • i Ci+l+J— p— 9) 

and similarly we can prove the formula is valid for j + I. Finally, we need to symmetrise the expression with respect 
to permutations of arguments C and C- o 

In the same manner, we can work out how a polynomial in u, v and their derivatives changes under the transformation 
T~^. Here we only give the formula. 
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Proposition 2. Under the transformation T , an arbitrary polynomial u^v-'bij{S,i, . . . . . . , Q) takes the Jo 



U V' J . . . , -^j., • . • , y 2*+^ 

p=0 5=0 

{{^iji^it ■ ■ ■ J Cp: Cli ■ ■ ■ I Ci-p: Cp+li • ■ • 7 Cp+gi Ci— p+li ■ ■ ■ i Ci+j— p— g)(Cl ' ' ' Ci+j—p-q) )c)C' 

where Cf are binomial coefficients. 

These two propositions immediately lead to the following result. 
Proposition 3. For any integer s > 1, if 

s s 



then 



■,s-l) 

i=0 1=0 



s min{i,l} 

h,s-ii^,, . . . , 6, Ci, . • . , Cs-i) = E E (19) 

i— p=max{0,l~s+i} 

{{bi,s-i{^l, • ■ • , ^p, Cli ■ • • ; Ci-p, Cp+1) ■ ■ ■ 1 6i Ci-p+1: ■ • • ) Cs-i)('?P+l ' ' ' OCi-p+1 • ' ' Cs-i)"^^ '')c)c 



and 



s min{i,j} 

Ml, . . . , c„ ci, . . . , c.-,) = ^ E E cf cn(-i)^-'^-^+^ (20) 

■ ■ • ,CpiCii • • ■ ,Ci-p,^p+i, ■ ■ • ,^j,Ci-p+i, ■ ■ ■ ,Cs-j){Ci ■ ■ -Cs-j) " ^'^'')i)<:, 

where Cf are binomial coefficients defined in Proposition^^ 

Proof. Applying formulae p7|) in Proposition [T] to every term in the sum and collecting the coefficients at every 
monomial u^v^ , we obtain 

s s i s — i 

J2 u^v'-%^s-Mi, . . . , ci, . . . , Cs-.) = E E E u'^^''i'-'^-''cfc!_,{-ir-^-'^ 

i=0 1=0 p=0 q=0 

{{bi,s—i{^l, ■ ■ ■ , £,p, Cl, ■ ■ ■ , Ci—p, Cp+lj ■ ■ ■ , £,p+q, Ci—p+1, ■ ■ ■ , Cs—p—q){£,p+l ' ' ' Cp+ijCi— p+1 ' ' ' Cs—p—q) '^^ ){)c 
s s min{i, 1} 

2—0 /=0 p=max{0,l — s+i} 

{{bi,s-ii^i, ■ ■ ■ ,^p,Ci, ■ ■ ■ , Q~p, Cp+ii ■ ■ ■ , ^i,Ci-p+i, ■ ■ ■ , Cs-i)(?p+i " ■ ■ 6Ci-p+i • ■ ■ ^)i)c, 

which leads to formula ([1^ . Similarly, using formulae (fTH]) in Proposition [H we can prove formula ^U^ . o 

Formula (I19p in Proposition [3] tells us how a polynomial of degree s changes under transformation T. The following 
result is an immediate corollary. 

Corollary 1. In variables u and v, system takes the form 

(4i+---+Ci+Ci+---+Cs-i) 2 (21) 

- _ 3 1 -i-s-i a,,,-,(gi,...,C.,Ci,...,C.-i) ' ' 

where ai^s-i{S,i, ■ ■ ■ ,£,i,C,i, ■ ■ ■ , Cs~i) defined in terms of aj,s-j{(,i, . . . , ^j, Cii • ■ ■ j Qs-j), j = ■ ■ ■ , s as in formula 
j'iP)) m Proposition\^ 



Proof. We perform transformation T to system \12\ : the linear part becomes diagonal and we simply write the 
system as 



Ut = U^r^ ' + ^^„+l-. (Ss>2 Ei=0 "'f'" 'a^,s-^ (6 , ■ • ■ , C^, Cl , • ■ • , C-*)) 

«t = -wC^' ;;^TT7 (Es>2E'=o'"*'^''"*a^'5~i('?i'---'C»'Ci,---,Cs-j) 



(22) 



Applying Proposition O wc obtain system (PT|) as stated. o 
For the future references we explicitly write down the relation between the quadratic terms of (fT^ and of (PTjl : 



02,0(6,6) 



«2,o(ei,6) + ao,2(6,6)(a6)"-*'^~'' + ^^01,1(^1,6)^2 " + 7T«i.i(6,Ci)C 



ai,i(a,6) = 2a2,o(6,6)-2ao,2(a,6)(Ci6)"+^"" + «i,i(6,6)^r" '"«i,i(a,6)^2"'"^ 
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00,2(6,6) = 02,0(6,6) + ao,2(6,6)(66)"+^"''-2«i,i(6, 6)6 

We now express the symmetries ([T^ and (IT5|) of system ([T^ in variables m, -0. 
Let us start with type I symmetries (flU) : 

s 

Mr = <r + .9, with 5r = ^^u*'y"-M,,,_,(6,...,6,Ci, 

s>2 2=0 

Due to the compatibility condition (jlip . its second component is of the form 

=<r + ^(A(5)), 

where the action of the operator Dt is defined in the symbolic representation as 



-2«ia(6,6)6^ 



, Cs — j)- 



(23) 
(24) 
(25) 



Dtig) = g*,n«) + ( + E E "'«'''o»,.-K6, . . . ,6, Ci, • • • , 6-^) 

s>2 i=0 



In variables w and v, a type I symmetry takes the form 



2i7"^5 



rf+2 + Dt 



(5) 
(5) 



(26) 



Here g stands for the transformed g under transformation T, whose exact expression can be obtained by applying 
Proposition [31 Notation Dt denotes the operator of total derivation with respect to t due to diagonal system (PT|) . 
Hence 

Dtig) = g*A{ut) + g*,v{vt). 



For type II symmetries (|15p : 

its second component is of the form 



with g ^^^u'v" Mi,s_i(6,...,6,Ci,---,6-i), 



s>2 1=0 



^ ^ «'t;^-*a,,,_,(6, . . . , 6, Ci, • ■ • , 6-0 + r?-''A(ff). 



, s>2 i=0 



Using the same notations as described for type I symmetries, in variables ii and v a type II symmetry takes the form 



7i-\-m-\-f- 



Dt 
-Dt 



g + iry™+'-"- 

n „ l„m+r—n 
y 2 ' 



' (Es>2 Ej=o "'"^^ *aj,5_,(6, • • ■ ,6,6,--- ,6-0 
Es>2Ej=o"*'^^~*^^'*-»(6,--- ,6,6,- •■ ,6-j)) 



(27) 



3.2 Approximate symmetries and integrability conditions. 

We derive the compatibility conditions of (fTi]) or (fT5|) and p2p via the compatibility conditions of their transformed 
forms (j26p or (|27p and (j2ip . For simplicity, we rewrite them in the following compact forms correspondingly: 





= ur!(6) H 


-el 


Vr 


= i;17(6)4 


-e2 


Ur 


- "f^(6) H 


-el 


Vr 


= -«f^(6) 


+ e2 


Ut 


= uw(6) + 




Vt 


= -WW (6) 





= E.>2 ELo '<.-2(6,---,6,6,---,6-0, J = 1,2 



p = Es>2 E»=o "'^^ >i,s-i(6, • ■ • ,6,6, ■ • ■ ,6-i)- 



(28) 
(29) 
(30) 
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Now we introduce a useful notation 

i s—i i s—i 

Gr,s-. ici,C2; Cl, ■ • ■ , C., Cl, • ■ • , Cs~^) = cMT. ^' + E Cfc) - ^ C^te) ^iCk). (31) 



1=1 k=l 1=1 k=l 



Proposition 4. System i28\) or system i29\) is compatible with i30\} of degree 2 if and only if the quadratic terms in 
or system i fl^j) are related to the quadratic terms in i30\) as follows: 



e2.ol4i'S2j — T^nm r-j — ^P2,ol?i,?2j, e2ol?i,?2j — / i 1 c — rT-P2,ol4i, ?2 j, 

t'2,oU7-i;?i,42j yo,2l~J^' ^;4i,?2j 

ei,i(4i,Ci) = T^HTTi 1 c /- x Pia(a.Ci), 61,1(6, Ci) = --;^:j-^—r-r^—p-:Pi,i{Ki,Qi), 

fe'l.iU) ~i; 4i, 4ij t/i,i(-i, i;Ci74ij 

1 ^ N _ ^"2(1: c; Ci,C2) //■ /- N 2 //■ ^ N _ g"o(c, i;Ci,C2) , , 

eO,2lU, S,2j — T^Hm , , > nPO, 211,1, (,2j, eg 2l<;i, 42; — ^/-w / 11/- /■ nP0,2(,(,1, <,2j, 

t'o 2U, 42j ^2 ol~^' -l^' "^1' S2j 



where c — 1 for the terms in \2^) and c = — 1 for the terms in 
Proof. Recall that for any polynomial pair 

F = M*t;''"*ai,s-i(6,-,Ci,Ci,--- ,Cs-i) and Q ^ u^v''^^hj^r~]{£.i,-,£.jXi, - ■ ■ Xr-j), 
where s>2,r>2,Q<i<s and < j < r are integers, we have [TT] 

The compatibility conditions of (^5]) and ([50]) up to degree 2 read as 

-«^(Ci) y'l,ELo"^«'-^2-.(a,---,e.Ci,---,C2-.) 



(32) 



Using formula (|32p and collecting the terms at every monomial for i € {0, 1, 2}, we obtain the formula for the 

quadratic terms in ()28p with c ~ 1. Similarly, we can derive the quadratic terms in ()29p as in the statement. o 

In fact, the results of Proposition |4] can be generated to the terms of higher degree i + j > 2, namely, one can derive 
that 

1 /t t /■ /• ^ (1, 61, ■ • ■ , Ci, • ■ • , Cj) /■ \ , tI roo\ 
eij(4i,---,4i,Ci,-- - ,0) = _,. . ■;rTP^,j(il,■■■,i^,Cl,■■■ ,Q) + Lij (33) 

2 f /• /• ^ 5j^i(c, 1; Cl, ■ ■ • , 0, Ci, ■ ■ ■ , 4i) ,^ ^ . /• ^ I r2 roA\ 

-l-, -"-I SI, ' ' ' ,S.j,Sl,''' 

where c = 1 for the terms in ([28|) and c = —1 for the terms in (|29p . Expressions Lj j and ^ depend only on the 
terms of degree lower than i + j in system pop . We will not write out them explicitly. 

Now we are ready to derive the compatibility conditions of (fn| or (|15p and Let us introduce the following 

notations: 

s s 

7vi™)(a,...,es) - (E^^r-E^r (35) 
1=1 1=1 

^1?(6,...,6) = (6 + ---+6)'"+^-| VC""- . (36) 



(E^r- 


- E C 




l=i+l 




E 



Mi:r^(a,...,e.) = (6 + ... + 6)" EC^^- E erM -EC^"^^+ E C^"^^- (37) 

\;=i ;=i+i / 1=1 i=i+i 
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Theorem 1. Expression Jj^p and expression il5]) are approximate symmetries of degree 2 of system ilS^) if and only 
if functions ^2,0; ^1,1; ^0,2 given by formulae 



^2,0 (6, 6) 

^0,2(6,6) 



4e 



2A2,o(6,6) + 2Ao,2(a,6) + ^i,i(Ci,6) + ^i,i(6,a 

2^2,0 (6 , 6 ) + 2io,2 (Ci , 6 ) - ii, 1 (fi , 6) - ^1, 1 (6 , 6 ) 



1 



are polynomials in ^1,^27 where ^2,0j ^0,2 are determined as follows: 
1. For the quadratic terms in \14^, 



(38) 
(39) 

(40) 



^2,0 (a, 6) 



(6+6)'w(™)(a,e2), 



'S'2"2'(^l'^2) 



02,0(6,6), ^1,1(6,6) 



(?i+6)'-A^i'"^(6,6). 



ai,i(6,6), 



^ . . (ei+6)''iV2^'"^(6,6) ^ 

^0,2(6,6) ^(„),. . , ao,2(4i,6)- 

'-'2,2 l?l,t2j 



2. For the quadratic terms in \15]) . 



(41) 



{ii + ^2YM!f^\ii,£,2) ^ (6 +6)'' ^27^6, 6) 

^2,0(6,6) = TTT— ^ — ^ 02,0(6,6), ^1,1(6,6) = TTT— ^ — ^ aia(6,6) 



^0,2(6,6) 



'5'2"2'(6,6) 

(6+6)''M^:^^(6,6) 

'5'2r2 (6,6) 



4^(6,6) 



00,2(6,6)- 



(42) 



Here 02.0, ai,i, 00,2 are given in terms 0/02,0,01,1,00,2 by i2^) . {24^ and \2S(I . 



Proof. Expression ([H]) and (fTK|) are approximate symmetries of degree 2 of (fT^ if and only if they are compatible 
with (fT^ up to degree 2. Therefore the diagonal forms of these symmetries ()26p and (f??!) must be compatible with 
the diagonal form of our system (PTjl . The latter is of the form of ([501 with Ci;(a;) = a;"+2, while and ([?7|l are of 



the form (US]) with Vl{x) = x"' and (HH) with f](a 



respectively. From Proposition [31 we know the relation 



of ^4.2.0, ^1,1, ^0.2 and j42.o, ^1.1, ^0.2 as well as the relation of 02,0, oi i, ao. 2 and 02,0, oi^i, 00.2- Now we need to show 
that vl2,o, ^1,1, ^0,2 are determined as stated. 

Consider first the case of type I symmetries ([H|) . Applying Proposition [D in the case of system ([50]) and ([^5[l with 
u>{x) = a;"+2 and Vi[x) = a;™, we obtain (for the coefhcient at v^v^): 



4,0(6,6) 



^^^'(6,6) 



(6+6)"+^-6^-C 



— P2,o(6,6) 



Comparing system ([^ to system ([^ and system ([5n[) to system ([^ . we have 



e^,o(6,6) 



(6+6)"+^+er^+er^ 



2(6+6)"+5 

Substituting these into the previous formula, we get 



^2,0(6,6) and P2,o(6,6) 



02,0(6,6) 

2(6+6)"+5-' 



. (6+6r^vf%^ 

^2,0ia,6j - ^ O2,o(,4l,42j, 

'5'2,2 (6,6) 

the first formula of ([^T|l . The other two relations can be obtained in the same way. 

Consider now the case of type II symmetries (|15p . Using Proposition [D in the case of systems pop and 
Lo{x) = and Vl{x) ~ nifi+'m+2 ^fJQ g^d (for the coefficient at v?v^): 

1 rc c \ (6+6) ^-?i -?2 c ^ 

62,0(6,6) = 1 — 1 — 1 P2,o(6,6)- 

(6+6)"+^-^r'-e2 ' 



with 



10 



Comparing system ([^ to system (P7)) , we know 



This leads to 

^2,0(^1,6) = "2,0(6, 6), 

"^2,2 l';i,S2j 

the first formula of (jH]). The other two formula can be obtained similarly. o 

The study of symmetries of ([T2|) is the study of the divisibility of certain special functions. If we introduce variables 
xi,X2 instead of by relations 

then functions 02,0(3^1: 2^2), ai,i(a:f , a;f ), aQ^2{x\,x\), S {x\,x'^) &ud M^\x\,X2) hecorae polynomials ina;i,X2(see 
formulae ^ and ^ ). ' 

Corollary 2. // and \15\) are approximate symmetries of of degree 2 then functions A2fl{xi, x\), Ai^i{x\, x\), 
and Aq^2{x\tx\) defined by CTl) and (OH) are polynomials in xi,X2- 



Proof. From Theorem [U if pi)) and are approximate symmetries of of degree 2 then 

^2,o(a;i, a;|), >li_i(a:^, a;|), and Aq, 2(2^1, 2^2) ^^'^ polynomials in cci, 0:2 (in fact these functions are polynomials in x\,x\). 
We know A2fi{x\, x\), Ai^i{x\,x\), Aq^2{x\^ x\) can be expressed as polynomials in ^2,0, ^14, ^0,2 with coefficients be- 
ing polynomials in a;i,a;2 according to Proposition [3l Therefore j42.o, ^1,1, and ylo.2 are polynomials in xi,X2. o 

We are able to obtain explicit recursive relations for determining the higher degree terms in (|14p and (|15p . Similar 
as formula (|33p and ((Ml) , we present the strictures of these expressions without the full details. From the proof of 
Theorem [TJ we only need to show how functions Ai_s-i are related to functions Oi^r-i, where r < s and these can be 
obtained by using formula (|33p and (p4|) and comparing the notations in (|28l) . (|29|) and pO|l with the ones in ([26]) . 
(EH) and (EB- 



Suppose that (HU or (fTK)) is an approximate symmetry of of degree p > 3. Then for any 3 < s < p and any 
i = 0, . . . , s we have 



Type I : . . . ,6) = tf^^^"^^^'' " " ' ■ ■ '^^^ + (^3) 



Type II: ^,^^-1(6, • ■ • , 6) = ■ ^ a^,s-^{^l, ■ ■ ■ ,(,s) + L'. ' , (44) 



where functions Ai_s-i and Oi^s-i are given in terms of Ai^s~i and Oi^s-i according to Proposition [31 and Lj-*-* 
depend only on lower degree terms a^j, i + j < s, and A^i'"\ S'^"^ and M^^™^ are defined as ([SS]) . ([36]) and ([37| . 

Again if we introduce variables 

1 1 

Xl — T ■ ■ ■ T X s — , 

functions ai^s-i{x\^ . . . ,x1), Ai^s-i{x\, . . . ,xl), S*^"' (xl, . . . , x^), and M^^™' (x^, . . . , a;^) all become polynomials in vari- 
ables Xl, . . . ,Xs. 

4 Global classification result. 

In this section, we state and prove the global classification result of system ^ in the sense that we list out all the 
integrable equations of arbitrary order. We begin with an crucial theorem leading to global classification. In order to 
prove the theorem, we prove the irreducibility of a family of special polynomials. 

Proposition 5. For any i — 0, . . . , s, polynomials 

s'-h^i ■■■,xi) = {xi+---+ xir^^ - +■■■+ x^+i - — -r^-'Y 

are irreducible polynomials in Xi, . . . , Xs over C when s > 3 and n > 1. 
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Proof. It suffices to prove that the polynomial si^} {xf, . . . , a;^), of. is irreducible over C since all the rest can 

be obtained from it by taking negation of s — i arguments Xi+i, . . . Xs- 

Let us first suppose that s > 4. If si"} is reducible then the projective hyperspace given by si^} — consists of two 
components. These components intersect in infinite number of points, which should be singularities of the hyperspace. 

(n) 

These singular points can be computed by setting the derivatives of Ss.J with respect to Xj, j = 1, . . . , s equal to zero. 
We obtain: 

^^^"^ - 2(2n + l)x, [xf +1 + ■ • • + xf+i) - (x? + ■ • • + = 0, J = 1, . . . , s. 



dxj 

So the singular points are either the solutions 

xl + ... + xl=0, +1 + ■ • ■ + xf+i = (45) 

or the solutions of 

xf-i = = . . . = xf-i = X, {xj + ■ ■ ■ + xlf"-' = x\ (46) 

The singular points from ([l5|) are kink points, which do not contribute to factorization. From it follows that 

{xlf"-^ = (x2)2"-i = . . . = (a;2)2«-i = {:,l + ... + xlf'^'K (47) 

Hence the coordinates differ by {An — 2)-th root of unity. Thus we get finitely many singular points. Thus the 
assumption that si"'} is reducible is false. 

When s — i, system (|Tf| writes as 

{x\f--^ = (x^)"""^ = {xir^' = {xl+xl+xir-\ 

By taking X3 = 1, we have that x\^ x\ and w are 2n — 1-st roots of unity such that w = x\-\- x\ + 1. Note that four 
complex numbers of the same absolute value can only add up to zero if they form the sides of a parallelogram with 
equal sides. Therefore either w — \ oy x^ ~ —1, i = 1,2. If w = 1, then x\ = —x\ leading to (xf)^""^ = — (a;|)^"~^, 
contradicting to {x\Y"-~^ = (a;!)^"^^- If a^f — —1, we have (a;f )^"^^ — —1, contradicting to {x\Y"~^ = (xg)^"^^ = 1. 
Hence, s'^l is irreducible. o 

The irreducibility of this family of polynomials leads to the following useful result in testing integrability. 

Corollary 3. If there are no quadratic terms in homogeneous system il^) with W{u) > 0, i.e. 

02,0(6,6) = 014(6, Ci) = 00,2(6,(2) = 0, 
and aij 7^ for some i + j > 2, then system il^) does not possess higher symmetries of any order. 

Proof. Without loss of generality, we assume that some of the terms in (|12p of degree s > 3 are not equal zero and 
all terms of degree less then s are equal zero. Suppose that such system possesses a type I symmetry (fTl|) of order 
m > 1. Then it is easy to see from formulae (j4T|) and (|43|l that yli.j(6, • ■ ■ , 6, 6, • • ■ , 6) = 0, i + j < s and 

A (r'^ 2n (^1 -I +x^)''ivi \xl,...,x'j,) 2 2^ ■ r. 

The left hand side of the above formula must be polynomial in xi, . . . , Xs. From Proposition[5]we know that polynomials 
5'|"^(xf , . . . , Xg), i = 0, . . . , s are irreducible polynomials over C and these polynomials do not divide 

ixl + --- + xlrNi"^Hxl...,xl) 

for any m > 1 and r > 0. Therefore S*^"^ must divide ai^s-ii^i, . . . , x'^), what is impossible since 

deg(a,,,_,(x2, . . . , x2)) = 4n + 2 - 2(s - l)W{u) < deg(^("^)(x2, . . . , x^)) = 4n + 2, z = 0, . . . , s. 

Thus system ([T^ does not possess any type I symmetry. The consideration in the case of type II symmetries is similar, 
o 

We are now ready to prove the following important theorem in classification of integrable homogeneous systems of the 
form H! 
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Theorem 2. Consider homogeneous system with W{u) > 0. Assume that it possesses an m-th order type I or 
type II symmetry of the form Jj^p or U5\). Suppose there is another system of the same weight and of the same form 

\v,^ + Y.u>2 Eto . . . , 6, Ci, ■ • ■ , a-^), ^^^^ 

whose quadratic terms equal to those of ilS^) . that is, bi,2~i{x,y) — ai^2-i{x,y), i = 0,1,2. Then if system 
possesses an m-th order type I or type II symmetry, then equation |^<$[ j and are equal and sharing the same type 
I or type II symmetry of order m. 



Proof. We will prove the statement by induction. Consider the case of symmetries of type I. Let us suppose that 
equation (j48p possesses a type I symmetry of the form 

s 

Ur = UC + E E ^^'^-^ (Cl , . . . , ^„ Cl , ■ • ■ , C.-, ) (49) 

s>2 j=0 

The quadratic terms of this symmetry coincide with those of ([Ti|) Bi^2-i{x,y) — Ai,2-i{x,y), i = 0,1,2 because the 
quadratic terms of and coincide. Let us prove that the terms of degree 3 in ^T7\i and (|15| coincide. According 
to formula P5)) we have 

A /2 2 2\ ("^1 I 2 i ^ 2 ) "^^^3 ( 1 ' 2 ' *^ s)'^ /2 2 2\ 7"(^) 

^i,3-i(a;i,X2,X3) = („N 2 2 2\ aj.3-i (a^i , , ) + Lj , 2 = 0,1,2,3 

'-'3,1 (^li 2^3) 

for the coefficients of the symmetry (ITi|) of ((T^ and 

B- , fT^ t2 - (^1 +^2 +3:^3) '^^3 \xi^X2iXl) ^ J 2 2 ^2x , f (3) „■ _ Q 1 2 3 

iii,3-i\Xi,X2,X^) — 2 2 2\ <}i^3-i{Xi,X2,X^) + 1^^ , 2 — U,i,Z,J 

for the coefficients of the symmetry ([49]) of (|48l) . 

Notice that Lf^ = \ i = 0, 1, 2, 3 since they only depend on the linear and quadratic terms of equations (|12p , (j48p 
and their symmetries (|14p and (HU, which are equal by the assumption. Therefore 

Bi^3~^i{Xi, X2, X^) — Ai^3^i(Xi, x\, x\) — 



\X\ X'2 yW^'^'ixl.xl^xl) / 2 2 2x ~ ! 2 2 2^\ 

— c(n) f 2 2 2\ \^i,i-i\Xl-,X2,Xj^) — ai^'i-i(Xi,X2,X3} \ . 

'-'3,1 V-^l' •'^2^ ^3) 

Its left hand side must be a polynomial in xi,X2,X3. In proposition we proved that polynomials s''^^^ {x\,x'^,x'^) 



-"3,1 \-^lT^2i -^d.) 

X2 

n > 1, TO > 1 and i — 0,1,2,3. Therefore S^"'^ {xf, X2, x^) must divide bi^3-i{xl, x^, x^) — ai^3_i(a;f , Xj, But this is 
impossible since, for i e {0, 1, 2, 3}, 

deg{k^3-,{xlxlxl) ~ a,^3^,{xl,xlxl)) = 2{2n + 1 - 2W{u) - r) < deg(5^"^ {xj^xl xj)) =4n + 2 
due to the assumption W(u) > and r > 0. Hence 

bi,3-i{Xi, X2, X3) — ai^3^i{Xi, X2, X3) = 0, Bi_3^i{Xi, X2, X^) — Ai_3^i{Xi, X2, X3) = 0, 
that is, fei,3-i(Cl,6:6) = ai,3-^i{£,l,£.2,^3), ^i,3-j(Cl, 6, 6) = A,3-i{£,l,^2,£.3)- 

Suppose we have proved that the terms of degree s — 1 coincide. For terms of degree s, we have, according to formula 

(BSD, 

; .2 2, ixl + --- + xlYNi"'\^u---,is) . ,2 2,^r(s) 

^i,s — i \Xi-, ■ ■ ■ ,Xgj — j.^^ ai^s — i\Xi,---,Xgj-\-Ln , I — U,...,S 



are irreducible polynomials over C. It is easy to see that it does not divide (x\ + X2 + x'^Y N^\x\,X2,X3) for any 



for the coefficients of the symmetry ([T4|) of (fT2|) and 

B -ix^ 2. {A + --- + xlYNi-^\i,,...,^,) ^ 

^J,s-jl,a;i, • ■ ■ , a-s j — TTT- ^ Oi^s^^\Xi, . . . ,Xg) + 1^,- , z — u, , 
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for the coefficients of the symmetry of (pS)) . Terms and Vf'^ depend only on coefficients of degree less than 
s and therefore they are equal i - — % z = 0, . . . , s. Hence 

f> (2 2\ _ A ( '2 {x\ + ■ ■ ■ X^Y Ns (^1 , . . . , ^s) /r I 2 _ ~ ('2 2\N 

i^i.s — i v^lT ■ ■ ■ T^s) ^i.s~i ■ ■ ■ i ^s) — (r,\ \ ^i,s—i[Xi, ■ ■ ■ , ) tti^s — iiXi, ■ ■ ■ , Xg ) j 



Again the left hand side of the last formula must be a polynomial in Xi, . . . ,Xs and therefore Sl"^/{xl, . . . ,x'^) must 
divide 5^^3-1(2^1, ■ • ■ ix'^) — di,s~i{x1, . . . This is impossible since 



deg(&i,,_,(x2, ...,xl)- a,^s-^ixl xD) = An + 2- 2{s - l)W{u) < deg(5^",^ (x?, . . . , xD) - 4n + 2, z = 0, . . . , s 

using the assumption W{u) > 0. Hence bi^s-i — o-i,s-i and Bi^s-i — for i = 0, . . . , s. Thus we proved the 

statement for type I symmetries. The proof in the case of symmetries of type H is analogous. o 

Finally, we state and prove the global classification result of homogeneous system ^ . 

Theorem 3. // a homogeneous system with Win) — w > Q possesses a hierarchy of infinitely many higher 
symmetries, then it is one 0/ the systems in the following list up to re-scaling u an, v — )■ /3f , t — > ^t, x — > 6x, where 
a, 13, 7, S are constant: 



Vt ^ U2+ iuvi + VUi — Su'^u 



2„. (50) 



vt = [D, + uf''{u) 71 = 1, 2, 3, ... . 



2n(„,\ „,2 ^ _ 1 O Q (51) 



Proof. We prove the theorem in symbolic representation. We know that system (jH]) in symbolic representation takes 
the form of (fT^ . In Corollary [3] we have proved that if (O possesses higher symmetries then the differential polynomial 
in ([9|) must contain quadratic terms and therefore w = W{u) > is either integer or half-integer. Indeed, 

T^ln+i-r+w = span{uiUj|i + j = 2n+ 1 - r - w, i,j e No}^ 

span{uiWj|i + j = n + i - w, i, j e Nq} ^ span{wjt)j|i + j = r ~ w, i,j e No}, 

where No = N1J{0}. Therefore if differential polynomial F in ([9|) contains quadratic terms then the weight of variable 
u is either integer or half-integer, i.e. either wGNotw = s— ^,s€N. Furthermore, from the Theorem [2] it follows 
that it suffice to classify integrable systems (0) up to quadratic terms. If system © or rather say possesses a 
hierarchy of infinitely many higher symmetries, it possesses infinitely many approximate symmetries of degree 2 of 
type I or type II. We first assume that these are of type I of the form 

«r, = <r + "'4'o(ei, 6) + uva[%^,(:^) + v^4%i,c2) + o{n^), j = i, 2, 3, . . . 

where 1 < mi < m2 < • ■ • < mj < TOj+i < ■ • • . From Theorem [T] and Corollary [51 it follows that 

aU) /■ 2 2\ / 2 I 2\r^^2 \-''1t^2)'- / 2 2\ Aij) / 2 2\ / 2 , 2\r-''2 \'^1t^2)'^ (2 2\ 

^2,0(2^1, 2^2) = (2^1 +2:2) — ^;^y— -— a2,o(a;i,a;2), AX;^{x^,x^) = [x^^^x^] ^ ax.i{x^,x^) 

'^2.2\!^\T^2) '-'2,1 ('''1 ' •'^2) 

A^^lixlxl) = (^^+^^)'^ \).7'^f ao.2(xj,x^) 

'-'2,2 (■^l! ^2) 

must be polynomials in 2:1, 2:2 for any mj, j — 1,2, . . .. Notice that polynomials 82"^ {x^, X2) and iVj™^"* (x^, Xj) can 
be factorized as 

'5'2,i (2^1, 2^2) = 2;]^X2>S'2 j (X]^, X2), N2 ' {Xi,X2) = X1X2N2 ' (X]^,X2), 

which leads to 

a(j) f 2 2\ I 2 , 2\r'-^2 1 1 2 ) '■ i 2 2\ A\3) I 2 2\ l 2 , 2\r^^2 1 J ''^ 2 ) - i 2 2\ 

^2,o(2;i,2;2) = (2^1 +X2) ^ 02,0(2:1, X2), ^l-;i(Xi,X2) = (2;i +X2) ^ Qi,i(Xi,X2) 



^0,21^1' 2^2^ — 12;^ + X2; —ao^2(Xi,X2 



^2,2 ('^1' ^2) *^2,1 (-^1' -^2) 

^2,2 (■^i I ■^2) 
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We can easily check gcd{{xl + x^Y , 82^^ {xf, xD) — gcd((a;? + x^Y , S2"'i {xf, xD) = 1. Moreover, from Lemma [T] 

in Appendix A it follows that if there exist infinitely many rrij such that Ai^2-i, i — 0,1,2 are polynomial, then 
rrij = 2, 3, 4, . . . and 

S2"'2ixi^xl)\a2fi{xl,xl), s!^''hxi^X2)\ai,iixl,xl), S^'^^ixf^xDlaoMxi^xl). (52) 
Counting the degrees of these polynomials we obtain that the division is possible only if 

deg(a,,2-») = 4n + 2 - 2r - 2w > deg(5^"^ (x?, xj)) = 4n - 2, i = 0, 1, 2 
implying r + w < 2. Since w > is either integer or half-integer, there are only four possible cases as follows: 

1. w = 2, r = 0; 

2. w = |, r = 0; 

3. w = l, r = 0, 1; 

4. w = i, r = 0,1. 

Let us compute the degrees of quadratic terms in the system we are considering. Applying formula (jl3p . we have 
deg(ai, 1(^1,^2)) = n - w + i, deg(ao,2(Ci, ^2)) ^ r - w, deg(a2,o(Cij C2)) ^ 2n + 1 - w - r. We know that aU these 
numbers should be positive integers, which enable us to determine the possible quadratic terms. 

Case[l] The only possible quadratic term in equation (fT2ll is 02.0. This leads to 

02,0(2;?, 2:2) = 00,2(2;?, 2:2) = 02,0(2;?, 2:2), ai,i(x?,X2) = 202,0(2;?, 2:2). 
From formula ([^^ we see that both 82^^ {x^, X2) must divide 02,0(2;?, 2:2). It is easy to check that 

gcd(5(:2^(x?,x2),^(;)(x?,x2)) = l. 

So we have S'2"2''(2;?,2;2)S'2"i (2;?, 2:2)102,0(2;?, 2:2). But this is impossible since 

deg(a2,o(2;?,a;^)) = 4n - 2 < deg(S'^"2^(x?, a;l)5^"i^ (x?, a;^)) = 8n - 4, for all n>0. 

Case [2] The only possible quadratic term in equation (I12p is oi.i. Notice that polynomial S'2"''(x?, X2), 02,0(2;?, x^) 
and 00,2(2;?, X2) are symmetric and polynomial ai,i(x?,x|) is antisymmetric with respect to xi,X2. Due to (|52p . we 
have 

02,0(2;?, X2) = -ao,2(x?,X2) = Ci(xi + X2)S^2'!2i^ 1^^2)1 ^l,! = C2(2;i - X2 )S';^"i (x? , X2 ) , (53) 

where Ci, C2 are arbitrary constants. Inverting formulae ([2111 and (|25p or using directly (|20p for ai,i, we find that 

01,1(2;?, X2) = 2^+1 (202,0(2;?, X2) - 2ao,2(x?,X2) + ai,i(x2,x?) - ai,i(x?,X2)) , (54) 
4X2 

which must be a polynomial in xi, X2. Substituting (|53p into (|54p . we see that this can only happen when ci — C2 — 
implying ai,i = 0. Hence our equation does not contain quadratic terms and therefore does not possess higher 
symmetries. 

Case [3] If r = 0, the only possible quadratic term in equation is 02,0- This leads to 

02,0(2;?, X2) = 00,2(2;?, X2) = O2,o(x?,x2), ai,i(x?,X2) = 202,0(2;?, X2). 
Similar to Case[Tl we obtain that S'2"2 (2;?, X2)>5'2"^ (x?, X2) must divide O2.o(x?,x2), which can only happen if 

deg(o2,o(x?,x2)) = 4n > deg(^("2^ (x?, x^)^'",) (x?, x^)) = 8n - 4, 
that is, 4n < 4 implying n = 1. Therefore, there exists a constant c such that 

02,0(2:?, x2) = c5^"2^(x?,x^)^^"^(x?,x2) = c (9xt + 9xt + Uxlxj) . 
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This leads to 

A^i^i^^l) = Ao,2{xl,xl) = c{3xl + 3x1 + 2x1x2) N^"''\ A[^\{xI,xI) = c(3x? + - 2xiX2)N^"''\ 
From formulae ^E^, and (gOl) we find that 

^2,o(-^i' ^2) = 3c(a;]^ + a;2)-/V2 ' (2^1:^2): ^1,1(2^11 2^2) = vIq 2(2;]^, X2) = 2c ■ 

X1X2 

It is easy to see that x\x\ does not divide N^'\x\,x\) for any integer rrij > 1, so Ao,2(a;i, Xj) is not a polynomial 
unless c = what implies 02,0 = 0. So system (1121) has no quadratic terms and therefore is not integrable. 

If r = 1, then we have oi^i = and therefore 

/ 2 2\ - / 2 2\ ^ 2 2\ I ^2 2\/ \2n+l 

a2,o(a;i,a;2) = ao,2(a;i, a;2) = a2,o(a;i, + ao^2(Xi, X2){xiX2) 
a'i,i{xl,xl) = 2a2fl{xl,xl) - 2aQ^2{xl,xl){xiX2f'^^^ 
It is easy to see that deg(oi,2-i(a;i, a;2)) = deg(S'2"^(x^,a;2)) = 4n - 2, i = 0, 1,2. It follows from ((S^ that 

a2,o{xl,xl) = aa.2{x\,xl) = ciS^2^l(xl,xl), ai,i(xl,xl) = 026'^"'' (x^, Xa), (55) 

where ci, C2 are constant. Inverting formulae P5)) . and (PS]) or using directly formula ((20)) for the quadratic terms, 
we find 

02,0(2;?, 0:2) = ^ {Aa2S){xl,xl) + ai,i{xl,xl) + ai^i{xl,xl)) , 

00,2(2:?, x^) = 2«-i 2»-i (402,0(2;?, a;^) - ai,i(x?,x^) - 01,1(2;^, a;?)) 



Hence 09,2(2;?, 2;2) = — 2ci. Without loss of generality we can put ci = ^, so 



Substituting (|55p into the above formulae we find that 00,2(2;?, 2;2) is a polynomial in x\,x^ if and only if C2 = 2ci. 

1 

2 ' 

2n 

ao,2(2;?,2;^) = -1, 02,o(:c?, 2;^) ^ -JZ^in+ii^ly-'i^l?''-' ■ (56) 

i=i 



Therefore 

^2"!o('^?' ^2) — ^o!2(^?i ^2) — 2 ^ ^2)-^2 '\xiix'^)i ^1^1(2;?, 2:2) = (2;? + a;2)Af2 ^''(2;?j2;2) 



and leads to 



^2'!o(^i'^2) — 2^x\+X2)N2 ^''(2;?, a;2), A[^\{xf,X2) — ^of2(^?; 2:2) — 0. 



Thus system (jl2p with quadratic terms (|56p and r = I possesses infinitely many type I approximate symmetries of 
degree 2 of any order mj > 1. 

Case [4l The only possible quadratic term in equation ([T^ is oi^i. If r = 0, then deg(ai^2--i(2;?, 2:3)) = 4n + 1, i = 
0, 1,2. Notice that polynomials S'2"''(x?, 2:3), 02,0(2;?, 2^2) and 00,2(2;?, 2:3) are symmetric while polynomial 01^1(2;?, 2:2) 
is antisymmetric with respect to xi, 2;2. Due to (|52p . we have 



02,0(2;?, 2:3) = -00,2(2;?, 2:3) = (ci(2;? +xl) +C2XiX2{xi +X2)) S^'^2i^l^^2)^ (57) 

ai^i{xl,xl) = {di{xl -xl) +d2XiX2{xi - X2)) Si^'^h^h^l)^ (58) 



where ci,C2,di,d2 are arbitrary constants. Substituting these expressions into ()54l) as in Case [21 wc find that 

4 2 

3C1, C2 = 3( 



oi,i(a;?, 0:2) is a non- vanishing polynomial in even powers of xi, X2 only when n = 1 and di = 2ci, (^2 = — |ci, C2 = |ci 



Without loss of generality we can choose ci = |. It follows that 

oi,i(a;?, 2:2) — II2;? + 92;2 

Therefore for Ai 2-i we obtain 



A^ilixlxj) - -Ai^lixlxj) = i^l{xl+xl) + x,X2{x,+X2)]Nt'\xlxl) 
A[]l{xl,xl) = (3(2;? - 2;^) - 2x1x2(2;! - X2)) Ni"''\xl,xl) 
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and hence 

0") 2 2 -^2"'^ (2^1,2^2) 2 2 0") 2 2 (j) 2 2 

^m(^1i^2) = ^ 2 (3a;2 +2X2) , ^2"|o(^l'^2) = ^0^2 (-^1' •''2) = 0. 



It is clear that A^"'|(a;f , a:|) is not polynomial since x\ does not divide -/V'2'"^'' (^^i; 2^2) for any rrij > 1. Thus equation 
()12p does not possess an infinite many type I approximate symmetries of degree 2. 

If r = 1, then deg(ai^2-i(a;f , Xj)) = 4n — 1, i = 0, 1,2. Similar as Case[21 we obtain the formula ([55]) . Inverting 
and I^S]) we find that 

ai,i(xi,X2) = 2!-! (2a2,o(a;?,a;2) - 2ao^2{xl,xl) +ai^i{xl,xl) - ai^i{xl,xl)) , (59) 
4a^2 

Substituting (j53p into ()59p . we find that ai_i(xf,a;2) is a non-zero polynomial in even powers of xi,X2 if and only 
C2 = 2ci and n ~ 1. Without loss of generality we choose ci = i, C2 = 1 and hence 



ai^i{xl, X2) — xf + 3x2- (60) 



For ■Al''l_^{xf, X2) we therefore obtain 



i^'i(x?,x^) = -A^^kxlxl) ^ l{x,+X2){xl+xl)Nt'\xlxl), 
A^Sixlxl) = {x,-X2){xl+xl)Nt^\xlxl) 

and hence 

4^(2^?, 2^2) = (x? +a:2)7v(™^)(a;?,a;^), ^2.0(2:?, 2;^) = Ao,2(a;^a:2) = 0. 

Thus we proved that system (I12p with the quadratic terms (|60p and n = — 1 possesses infinitely many type I 
approximate symmetries of degree 2 of any order rrij > 1. 

Summarizing the results in the above four cases, we proved that if system ([9|) possesses infinitely many type I ap- 
proximate symmetries of degree 2 of any order rrij > 1, then it is up re-scaling one of the systems in the following 
list: 



Ut = Vi 

Vt = U2 + uiv + ?>uvi + h, h = 0(7?.^) 

Ut = Vi 



2^ (61) 



Vt ^ U2n ~ + ^Y.'jll^2n+l'^3-l'^2n-j + h, h = o{'R^) ^^'^^ 

Let us now assume that system (jl2p possesses infinitely many approximate symmetries of degree 2 of type II of the 
form 

Ur, = + U24^^(a, 6) + Ut'^-lkCl, Cl) + f'^0.2(Cl, C2) + 0(7^'), J = 1, 2, 3, . . . 

and 1 < mi < m2 < ■ ■ ■ < rrij < mjj^i < ■ ■ ■ . From Theorem [T] and Corollary [21 it follows that 

'-'2,2 (2^1 ' 2^2) '^2,1 i-^i , 2:2) 

4(i)^^2 2n _ / 2 , 2Nr ^''-'2.2 ^^I'^zJ .2 2x 

'-'2,2 ("^i' ■^2) 

must be polynomials in xi,X2 for every rrij. Notice that polynomials M2^^'\x\,x\) can be factorized as 



M!^^\xW2)-xlxim\xlxl). 



which leads to 



m(™'3)(^2 ^2^ j('Amj),2 „2\ 

AVJl (^-^ ^-^^ _ (^-^ , r^Y^jS^_ril^n, Jr^ r'^] (r^ - (r^ . r'^Y {^1^^2) ,2 2n 

^2.'Q\^\i^2) ^ \X\+ X2) .(-^N 2 2 02,01,2^1, a;2 J, A^y^^X-y,X2) — \Xx+ X2) ^^^s ^ ^ ai,i(^Xi,X2j 

'-'2,2 l-^l ' 2^2 J '^2,1(2^112^2/ 



^0,2 12^1, 2:2; — I2;i-t-a;2j ^(„) ^ ^"0,21,2:1, a;2j 

'-'2,2 (,2^1 1 2^2/ 
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from Lemma [2] in Appendix A, if there exist infinitely many rrij such that the above expressions are polynomials in 
xi,X2, then rrij = 1, 2, 3, 4, . . . and 

S2'^2ixl,xl)\a2Mi^l,xl), s!£^{xl,xl)\ao^2{xl,xl) S!^"l{xl,xl)\ai^i{xl,xl). (63) 

We then repeat what we did for the case of type I symmetries and we obtain the similar conclusion: If system (|12p 
possesses infinitely many type II approximate symmetries of degree 2 of any order rrij > 1, then it is up re-scaling 
either (|6T|) or ([62ll . 

By now we have classified all systems (jl2p possessing infinite many approximate symmetries of degree 2. We know 
that system ((5(1)) possesses infinitely many higher symmetries of any order of form (fTH) and (fT5)l . Notice that the 
quadratic terms of (|50p and of the system (pTjl coincide. By Theorem^ if system (|FT|) possesses an exact symmetry of 
some order m for some h in its right hand side , it must coincide with the system (|50p . that is, h = —3u^ui. Similarly, 
quadratic terms of ([5T|) and of the system coincide and system ([51]) possesses higher symmetries of any order. 
Therefore, if system (|62p possesses an exact symmetry of some order for some h then it coincides with ([5T|) . o 

System (jSip is equivalent to a linear equation wu = d'^"^^w under the Cole-Hopf transformation u = (logw)^. This 
system possesses infinitely many type I symmetries 

u^^ = D^{D^ + u)"'-^u, m = 2,3, . . . 

and type II symmetries 

These symmetries correspond to the symmetries = Wm and = Wtm of Wu = W2n+i- 

5 Conclusion and Discussion 

This is the first paper dealing with the global classification of the systems with two components. The success of global 
classification relies on so-called implicit function theorem in 11]. The theorem implies that the existence of infinitely 
many approximate symmetries of low degree together with one symmetry leads to integrability. The proof of the 
theorem is straightforward and purely algebraic. However, to check the conditions of the theorem requires to prove 
the irreducibility of a family of polynomials. In this paper. Theorem [5] can be viewed as another version of it. The 
required irreducibility of polynomials S^^^ {x\, . . . , x^) defined by (j36p is proved by algebraic geometrical arguments in 
Proposition [5] 

The formal diagonalisation approach proposed in this paper can be applied to large classes of systems with the 
diagonalisable linear part. This type of equations is generic and it has not been treated by other classification 
methods. Our approach enables us to explicitly write down the necessary conditions for integrability. We are able to 
prove that these conditions are indeed sufficient. Resolution of these conditions leads us to the ultimate classification 
result. 

We have not applied this approach to the systems with nilpotent linear part yet. However, we do not expect any 
difficulties in this case. The final goal will be to combine these two classes and to obtain the global classification 
of homogeneous polynomial evolutionary systems of two components as it has been done for the scalar evolutionary 
equations [TT] . 

Appendix A 

Here we prove the lemmas required in the proof of Theorem [31 which is based on the theorem of Lech-Mahler and its 
straightforward corollary. 

Theorem 4 (Lech, Mahler). Let Ai, A2, . . . , An G C and ai, a2, . . . , a„ be non-zero complex numbers. Suppose that 
none of the ratios Ai / Aj with i ^ j is a root of unity. Then the equation 

aiA™ + 02^2"' + ■ • • + a„A™ - 

in the unknown integer m has finitely many solutions. 

Corollary 4. Let A, B,C £ C and a, b, c be non-zero complex numbers. Suppose that the equation 

aA"" + bB"" + cC"' = 
has infinitely many integers m as solution. Then A/B and A/C are both roots of unity. 
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We simplify our notations by writing S+ ^ S^^^^, S- = S"^"/, = M^"^'' and M„^. = M^"^^ 
Lemma 1. Consider an infinite sequence of polynomials 

^2™' (2^11 ^2) = (xf + xj)"'^ - xl"^' - xl"^' , 1 < mi < m2 < ■ ■ ■ < rrij < nij+i < ■ ■ ■ (64) 

// there exist infinite many nrij such that each N^^'^x^^x^) has nontrivial common divisor with the polynomial 
S'±(x?,xi) defined by (EB, Then ^ 2,S,A, ■ ■ ■ and gcd{S^{xl, xj), N^^\xl, xj), N^^\xl, xj), . . .} ^ xfxl. 

Proof. Since the polynomials S^ixf, x\) and N^'\x\^ x\) are homogeneous polynomials in xi, X2 it is convenient to 
consider them in afhne coordinate If they have nontrivial common divisor, these exist g as a common root, that 
is, 

S^{q) = (1 + g2)2n+l _ ± ^2„+l)2 ^ q 

^2"'^ (9) = (1 + 9^)™' - 1 - 9^"' = 0- 

which is required to be satisfied for infinitely many integers ruj > 1. Applying the Lech-Mahler theorem to it we find 
that this is possible if one of the listed cases holds: 

• (7 = (double root), ruj = 2, 3, 4, . . ., 

• q = ±«, nij — 1 mod 2, 

• q — e~ , e~3~, e~3~, e^, nij = 5 mod 6, 

• q — e^, e^, e^, (double roots), nij = 1 mod 6. 

The last three sets of roots do not satisfy the equation S^{q) = for any n > 0, while g = is a double roots of 
S^{q). Lemma is proved, o 

Lemma 2. Consider an infinite sequence of polynomials 

M± (x?,x2) = ixl+xl){xl-+'±xl"+')-xl-+'"^+' ~ i±xr+^"^+'), (65) 
< mi < 7712 < ■ ■ ■ < rUj < TOj_|_i < • • ■ 

If there exist infinite many rUj such that each M^.{x\,x\) (or M^.{x\,x\)) has nontrivial common divisor with the 
polynomial {x\,x\) (or {x\,x\) ) defined by \3b]) . Then mj — 1,2,3,4,... and 

gcA{S+ {xl , x2) , M+ {xlxl),M+{xlxl),...} = xjxj ; 
gcd{ 5- {xlxl),M^{xl,xl),M^{xlxl),...} ^xlxl 

Proof. Similar as in the proof of Lemma [l] we consider polynomials S^{xi,X2) and AI^.{x\,x'\) in affine coordinate 
1^. If S^{x\^x\) and M^.{x\,x'^) have nontrivial common divisor, these exist g as a common root, that is, 

S+{q) = (1 + g2)2"+l _ 1 _ qir^+2 _ ^^2n+l ^ q 

[q] = (1 + q^r^ (1 + 9^"+!) - {q^r^q^-+' -1 = 0, 
which is required to hold for infinitely many positive integers ruj . We apply the Lech-Mahler theorem to it and obtain: 

• q = (double root), ruj G Z+, 

• q = e~,e'^ and n — 1 mod 3, nij = mod 3 or n = 2 mod 3, nij = mod 6, 1 mod 6 or rt = 
mod 3, ruj = mod 6, 5 mod 6, 

• q — e~5~,e~3~ and n = 1 mod 3, rrij = mod 6 or n = 2 mod 3, rUj = mod 6, 4 mod 6 or n = 
mod 3, nij = mod 6, 2 mod 6. 

It is easy to check that the roots from the last three cases are not the roots of S~^{q) for any integer n > 1, while 9 = is 
a double root of S~^{q). Therefore, we have rrij = 1, 2, 3, . . . and gcd{S~^ {x^, x^), Mj^(a;f , a;^), M2*"(a;2, a;^), • • ■} = a;^. 
Similarly, we can prove the second half of the lemma, o 
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